Transport coefficients near chiral phase transition 
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We analyze the transport properties of relativistic fluid composed of constituent quarks at finite 
temperature and density. We focus on the shear and bulk viscosities and study their behavior near 
chiral phase transition. We model the constituent quark interactions through the Nambu-Jona 
Lasinio Lagrangian. The transport coefficients are calculated within kinetic theory under relaxation 
time approximation including in-medium modification of quasi-particles dispersion relations. We 
quantify the influence of the order of chiral phase transition and the critical end point on dissipativc 
phenomena in such a medium. Considering the changes of shear and bulk viscosities along the phase 
boundary we discuss their sensitivity to probe the existence of the critical end point. 

PACS numbers: 25.75.Nq, 24.85.+p, 12.39.Fc 



1. INTRODUCTION 

The shear r\ and bulk £ viscosities are parameters 
that quantify dissipativc processes in the hydrodynamic 
evolution of fluid. Furthermore, it has been argued 
that £ and jn are sensitive to phase transitions in a 
medium [J 0, i, |, i, % 0, @, ©, Q3. Indeed for cer- 
tain materials, e.g. helium, nitrogen or water, the shear 
viscosity to entropy ratio r\j s is known experimentally to 
exhibit a minimum at the phase transition [f|. On the 
other hand, the bulk viscosity C/s was argued to be large 
at the critical point 0, H, H, [lQ| • 

In QCD of particular interest are properties of trans- 
port coefficients near the critical line where deconfinc- 
mcnt and chiral phase transition sets in. This is because 
any change of the bulk and shear viscosities near T c modi- 
fies the hydrodynamic evolution of the QCD medium and 
influences phcnomcnological observables that character- 
ize its expansion dynamics. 

The recent Lattice Gauge Theory (LGT) calculations 
in pure SU(3) seem to be consistent with the expectation 
of decreasing r//s and increasing £/s toward the first or- 
der deconfinement phase transition approaching from the 
high temperature phase [ll|, G3 ■ The properties of the 
energy-momentum tensor in the vicinity of the second 
order phase transition and its importance for the anal- 
ysis of transport coefficients and their critical behavior 
have been also calculated and discussed within LGT for 
(3+l)-dimensional SU(2) gauge theory [l3|. All these re- 
sults show that the transport coefficients are of particular 
interest to quantify the properties of strongly interacting 
relativistic fluid and its phase transition [5(. 

In the theory of critical phenomena, describing dy- 
namical processes near critical point and their universal 
behavior, the transport coefficients in the second order 
phase transition can be quantified by a dynamical criti- 
cal exponents z [Til . Il5j . These exponents are common 
for all models belonging to the same dynamical universal- 
ity class. In general, the critical exponent z determines 
the critical slowing down of the system relaxation time 



t ~ £ z ru \T — T c \~ vz near the critical temperature T c 
with £ being correlation length and v its static critical 
exponent. It was argued that within the classification 
given by Halperin and Hohenberg the QCD critical end 
point (CEP) belongs to the universality class of the model 
H [lj| . Consequently, both shear and bulk viscosities are 
expected to diverge at the CEP as rj ~ and ( ~ £ z <-' 
with dynamical critical exponent z v ~ 1/19 and zq ~ 3 
respectively [lj|. 

The analysis of experimental data obtained in heavy 
ion collisions at RHIC showed that the evolution of the 
quark gluon plasma (QGP) is well described by nearly 
ideal hydrodynamics [17|. It has been even argued that 
the QGP is in fact the most nearly ideal fluid known, 
with shear viscosit y cl ose to the conjectured lower bound 
in any system 0,13- Clearly, in the evolution of QCD 
medium its transport properties are changing when ap- 
proaching the phase boundary. Thus, it is of importance 
to quantify the change of shear and bulk viscosities with 
thermal parameters near the phase transition expected 
in strongly interacting medium. 

In this paper we explore the transport properties of 
relativistic fluid composed of constituent quarks at finite 
temperature and density. The dynamics of medium con- 
stituents is described by the chirally invariant four-quark 
interactions within Nambu-Jona-Lasinio (NJL) model 
[20l . |2lj . The model correctly describes the universal 
critical behavior of physical observables with static criti- 
cal exponents that belongs to the QCD universality class 
with respect to chiral symmetry 

To calculate the change of the shear and bulk viscosi- 
ties with thermal parameters, we assume that the system 
appears near equilibrium and we apply the kinetic th eory 
under relaxation time approximation [22|, [H, HH US Hfif. 
The relaxation time of quarks is obtained from the ther- 
mally averaged total cross-section of elastic scattering in 
the dilute gas approximation 0]. 

In the NJL model under mean field dynamics the in- 
teractions of quarks lead to the quasi-particle description 
of thermodynamics with their masses M(T, fi) which are 
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temperature T- and density /i-dcpcndent [2l|, |27j ■ The 
dynamically generated mass acts as an order parameter 
for chiral phase transition. In our calculations of trans- 
port coefficients, following the method described in Ref. 
[24| . the constituent quark mass M(T,fi) is consistently 
built in. We quantify the influence of the order of the 
chiral phase transition and the critical end point on the 
shear and bulk viscosities. We show that the bulk vis- 
cosity is strongly increasing whereas the shear viscosity 
is decreasing when approaching the CEP from the side 
of chirally symmetric phase. These results are consistent 
with the expectation that the bulk viscosity could domi- 
nate dissipative hydrodynamics near CEP [2,11]. Consid- 
ering the changes of r\ and ( along the phase boundary 
we discuss their sensitivity to probe the existence of the 
critical end point. 

Within the kinetic theory and under linear response 
and relaxation time approximation used in this calcula- 
tion, there is no access to the critical dynamics near sec- 
ond order transition where the properties of shear and 
bulk viscosities are quantified by dynamical critical ex- 
ponents. In addition, in this formulation, the transport 
coefficients are finite with respect to static critical ex- 
ponents [24[. Thus, the behavior of r] and £ is entirely 
governed by the non-singular part of the partition func- 
tion. Consequently, the change of the shear and bulk vis- 
cosities with thermal parameters obtained in this paper 
may depend on specific model dynamics. Nevertheless, 
one expects that some of our results, e.g. that on the 
influence of the order of the chiral phase transition on r\ 
and £, could be relevant for qualitative understanding of 
transport properties in the QCD medium. 

The paper is organized as follows: In Section [2] we 
introduce the NJL model and its thermodynamics. In 
Section [3] we describe the methods used to calculate the 
shear and bulk viscosity. We also quantify in this section 
the properties of the shear and bulk viscosity near the 
phase transition and along the phase boundary. Finally, 
in Section [5] we summarize our results. 



2. THE TWO-FLAVOR NJL MODEL 

In order to quantify the transport properties of 
fluid composed of constituent quarks near the chiral 
phase transition we use the Nambu-Jona-Lasinio (NJL) 
model [2p| . This model exhibits basic properties ex- 
pected in QCD due to the chiral symmetry restoration. 
We start with the following Lagrangian for two quark 
flavors and three colors [211 



C 



G, 



(2.1) 



where to = diag(m tt , and fi = diag(/i u , fid) are the 
current quark masses and the quark chemical potentials 
respectively, whereas f are Pauli matrices. The coupling 
of four-fermion interactions GsA 2 = 2.44 and the three 
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FIG. 1: The phase diagram of the NJL model. The full and 
dashed lines show the first order and cross over transition re- 
spectively. The dot indicates the position of the CEP located 
at (T, fx) = (80.9, 329) MeV. 



momentum cut-off A = 587.9 MeV are chosen so as to 
reproduce the vacuum pion decay constant and chiral 
condensate with m u = = 5.6 MeV. 

In the mean field approximation and for the isospin 
symmetric system the thermodynamics of the NJL model 
is described by the potential: 



n{T,/j,;M)/V 



(M - to) 2 



12 



d 3 p 



[E(P) 



AGs " J (2tt) 3 
- rin(l - f(p, T, n)) - TIn(l - f{p, T, fi)] , (2.2) 



with M = to — 2Gs(ipip) being a dynamical quark mass, 
E(p) = V? 2- 



M 2 the quasi-particle energy and /, / 
-l 



1 + exp [(E(p ) =F /x)/T] J are the particle / antiparticle 

distribution functions. The dynamical quark mass M in 
Eq. (|2.2|) is obtained self-consistently from the stationar- 
ity condition dCl/dM = that leads to: 



M = to + 12G S 



d 3 p M 
(2tt) 3 ~E 



[1-/-/] 



(2.3) 



The thermodynamic potential of the NJL model ex- 
hibits a typical QCD-like phase diagram shown in Fig.[T] 
The cross-over transition is identified from the peak po- 
sition of the chiral susceptibility 



Xch 



d(tpip) 
dm 



(2.4) 



The CEP is obtained from the condition of the zero- 
curvature of an effective potential or equivalently from 
the divergence of the quark-number susceptibility. 

Our main goal is to quantify the transport properties 
of a medium composed of dynamical quarks with mass 
M(T, fj) near the chiral phase transition and along the 
boundary line shown in the Fig. ([T]). In this context we 
explore the shear 77 and bulk £ viscosities. 
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In the thermodynamic system 77 and £ are sensitive 
to thermal and dynamical properties of a medium. To 
explicitly calculate transport coefficients we use the ap- 
proach based on the relativistic kinetic theory. 



coefficients of the space-space component of the energy- 
momentum tensor away from equilibrium. 



3. TRANSPORT COEFFICIENTS 

In the relativistic kinetic theory the transport param- 
eters, the shear and the bulk viscosities, are defined as 



For the fluid composed of quasi-particles with dynam- 
ical mass M(T, fx) and under the relaxation time approx- 
imation,};]^ shear and bulk viscosities are obtained as 
follows m 



d 3 p p 1 
5T J (2tt) 3 E 2 
A f d 3 p M 2 
Tj (2^)3^ 



[r/(l-/)+f /(!-/)] 



(r/(l-/)-r/(l-/))( — 



(r/(l-/)+f /(!-/)) (j^ 



dP 



dP 
Ik 



E-T 



BE dE 



dT 



/' 



dP 
dn 



dE 



(3.1) 



(3.2) 



The derivatives of pressure P with respect to the net 
quark number density n or energy density e in Eqs. (|3.1j) 
and (|3.2[) can be expressed in terms of susceptibilities 
Xxy = d 2 P/dxdy as 



dP 
Ik 
dP 
dn 



uTxtt + (n/i 



?T)x^T - s^Xnn 



(3.3) 



with Cy being the specific heat calculated at a constant 
volume. The Cy is also expressed through susceptibilities 



C v =T 



ds 
dT 



T 



Xtt 



xIt 
Xfj.fi 



(3.4) 



To quantify the change of the above transport coeffi- 
cients in quasi-particle medium one still needs to spec- 
ify the relaxation time r. The r is calculated from the 
thermal averaged cross-section describing the total elastic 
scattering of medium constituents. In the NJL model the 
medium is composed from constituent quarks with two 
different flavors. Thus, in the dilute gas approximation 
the relaxation time r is obtained from [19, 281 



t 1 = 6 [n q (a 

+ nq (a. 



uu — >uu 



I" &ud->ud) 

+ &uu^dd + Vud^ud)] > (3-5) 



where n% are the quark/antiquark densities 

d 3 p 



(2tt) 3 
d 3 p 
(2^)3 



f(p; t, fi) , 
f(p; v) ■ 



(3.6) 




FIG. 2: The collision time in the NJL model at /i = 0. 



The in-medium cross section for quark- antiquark, 
quark-quark and antiquark-antiquark scattering pro- 
cesses were studied in detail in Rcf. d including 1/N C 
next to leading order corrections. For uu — > uu scatter- 
ing the differential cross section is obtained from the T 
matrix as 



d>0~uu — >uu 

dt 



(s,t;T, u) 



16tts(s- AM 2 ) 



\T U u^uu\ 2 (s,t;T,fx) , (3.7) 



where the average over the initial color and spin as well 
as the summation over final states were taken. We con- 
sider the rest frame of colliding quarks so that the differ- 
ential cross section (|3.7[) is expressed as the function of 
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the Mandelstam variables s and t. The integrated cross 
section is calculated from 



dt 



da U u->uu -4t(s + 1 - AM 2 ) 



dt (s - AM 2 ) 2 

x (1 - /(Vi/2 - M )) (1 - /(V?/2 + fi)) . (3.8) 

In the above equation we have taken into account that 
the transpo rt p rocess is dominated by the large angle 
scattering j28[. The terms in the brackets correspond 
to the Pauli blocking factors which account for a pos- 
sible occupation of particles in the final state. We also 
introduce the energy-averaged cross section 



' uu — > uu 



(2» 



dso, 



uu — >uu 



[a;T,n)P(a;T,n), (3.9) 



with P being the probability of yielding qq pair with the 
energy s 



P(s; T, n) = C^fs{s-AM 2 ) 

x /(VS/2 - n) f{yfs/2 + fi) v Te y(s) , (3.10) 

where the relative velocity between two particles in the 
initial state v re i is given by 



Viel(s) 



AM 2 



(3.11) 



The normalization constant C is fixed from 

J dsP(s;T,n) = 1. (3.12) 

The differential cross sections for different quark flavors 
are obtained through the isospin symmetry, charge con- 
jugation and crossing symmetry. The total cross section 
after the energy averaging reads 



*dd 



id — >ud * ®uu — >uu ~r @ud — >ud • 

(3.13) 

The dominant contribution to the total cross section 
(|3.13p comes from the s-channel due to the propagation of 
7r and cr modes. The mass of a drops toward the critical 
temperature T c and due to chiral symmetry restoration 
is degenerated with the pion mass at T = T c . On the 
other hand, for T > T c , both tt and a masses are increas- 
ing with temperature. Consequently, the cross section 
a(s = m^iT, fi) shows singularity at T c 3 # x . How- 
ever, after the integration over energy s this singularity 
is mostly washed out in the total cross section a q . As 
remnants of the above singularity a broad bump in the 
a q is seen arround T c [3[. Thus, the total cross section is 
non-singular even at the CEP. 



Fig. [2] shows the collision time calculated from Eq. 
(|3.ip at finite T and at /j = 0. There is a clear decrease 
of r when approaching the critical region from the chi- 
rally broken phase. One should realize that in spite of the 
phase transition the collision time is quite smooth around 
T c and experiences only a shallow minimum at the phase 
transition. This behavior is due to the total cross section 
which shows a similar broad structure around T c . The r 
is large at low and small at high temperatures since the 
quark density is increasing function of temperature 

At the finite fi the total cross-section is larger than 
that at fi = and the quark density increases while the 
anti-quark density decreases. Consequently, the collision 
time for quark r and for anti-quark f behaves as 



r(/0 > r(ji = 0) , t{h) < t((i = 0) 



(3.14) 



With the above results for the collision time applied 
in Eq. (|3.1[) we can calculate the transport coefficients 
in the quasi-quark medium and study their sensitivity to 
the chiral phase transition. 



3.1. Bulk and shear viscosities near phase 
transition 

The temperature and density dependence of the shear 
viscosity obtained in the NJL model under relaxation 
time approximation is shown in Fig. [3] For low den- 
sity the shear viscosity shows a shallow minimum near 
the pseudo-critical temperature. This behavior is similar 
as that observed in the collision time. The constituent 
quark mass decreases with temperature and the quark 
density is enhanced in the chirally symmetric phase which 
results in the rise of rj at high temperature. Quantita- 
tively the difference in i] at /i — and ii ^ comes mainly 
from the density effect. Consequently, at the CEP and 
near the first-order phase transition, this effect should be 
even more significant. Fig. [3]- (right) shows the change of 
shear viscosity when crossing the CEP and the boundary 
of the first order phase transition. As expected, the 77 
shows discontinuity at the first order transition. At the 
CEP this discontinuity disappears, however there is still 
an abrupt change of 77. One observes that the shear vis- 
cosity is decreasing toward the CEP and first order phase 
transition from the side of the symmetric phase. 

When discussing the transport properties of relativis- 
tic fluid a natural normalization of shear viscosity is the 
entropy density s. In this way the rj/s is dimension- 
less quantity that characterizes dissipation of energy in 
a medium. Fig. [4] shows the temperature dependence of 
the ratio rj/s at vanishing and at finite chemical potential. 
The entropy density rapidly increases with temperature, 



# x More precisely, the singularity appears at the pionic Mott tem- 
perature at which the pion mass becomes equal to the mass of 
its constituents, = 2M. In the chiral limit the Mott temper- 
ature coincides with the chiral phase transition temperature. 



* 2 The behavior of t at T < T c seen in Fig. [2] is specific for the 
NJL model where in the low-temperature phase the medium is 
composed only of heavy constituent quarks. 
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FIG. 4: The ratio of shear viscosity to entropy density around cross over (left-hand figure) and around /xqep and first order 
transition (right-hand figure). 



thus the rj/s is reduced in the high-temperature phase. 
A shallow minimum structure around T c for small fx is 
not preserved with increasing fi. From Fig. 2]- (right) one 
sees that rj/s is insensitive to the existence of the CEP. 
At higher [i the ratio shows a jump associated with the 
first order phase transition. 

From the general analytic expression for the shear and 
the bulk viscosities (see Eqs. (|3.2j) and (|3.ip ) one expects 
that £ should be more sensitive to the phase transition 
than r\ due to explicit contribution of derivative terms 
involving a dynamical quark mass. Fig. [5] shows the tem- 
perature dependence of £ for different (i. In the regime of 
cross-over transition the density dependence of £ is much 
weaker than that seen in the shear viscosity. This is in- 
deed due to dM/dT and dM/dji terms that contribute 
to the bulk and are absent in the shear viscosity. The 
temperature dependence of £ at finite \x is similar to that 



at n = as seen in Fig. [5]- (left). With the increasing 
chemical potential, the quark density effects start to be 
visible. The growth of £ toward T c in Fig. [SJ- (right) is 
due to an increase of quark number density. At the CEP 
the bulk viscosity is finite in spite of the fact that all sus- 
ceptibilities contributing to C are diverging (24|. Thus, 
at the CEP there is a precise cancelation of singularities 
from different terms in Eq. I|3.2|) . Nevertheless, one ob- 
serves a rapid downward change approaching the CEP 
from the low-temperature side. This reflects the change 
of the dynamical quark mass with temperature. 

The ratio £/s, of the bulk viscosity to the entropy den- 
sity, keeps basically a structure of (/A 3 . The suppressed 
entropy density due to the heavy constituent quark mass 
in low-temperature phase results in growth of the ratio 
shown in Fig. [5] The rate of change in Q/s with tem- 
perature at CEP is higher than that of rj/s. However, 
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T/T c T/T c 

FIG. 6: The ratio of bulk viscosity to entropy density around the cross over (left-hand figure) and around hcep and first order 
transition (right-hand figure). 



the absolute value of (/s is smaller than rj/s at CEP. 
This result is valid under relaxation time approximation 
where r is always finite. However, when going beyond 
this approximation and including dynamical scaling [l5| , 
the bulk viscosity should be singular at the CEP. 

In the derivation of the bulk viscosity we have explic- 
itly included in Eq. (|3.2|) the modification of particle dis- 
persion relations due to a dynamically generated quasi- 
particle mass which is T and \i dependent. In Fig. [7] 
we quantify the influence of derivatives of a dynamical 
mass on the bulk viscosity. The full-line in this figure 
was obtained from Eq. (|3.2p whereas the dashed-line was 
calculated by eliminating the dM/dT and dM/d/j, terms 
from Eq. (|3.2[) . The behavior of ( is largely governed by 
the dynamical quark mass M(T,fi) multiplied in front 
of (|3.2p . Unlike the expectation that the bulk viscosity 
should be sensitive to the phase transition, there is no 



much difference in the above two cases. The enhance- 
ment of dM/dT around T c is weakened by the pre-factor 
M 2 in Eq. (|3.2j) . In addition, the contributions of differ- 
ent susceptibilities in Eq. (|3.2p which are diverging at the 
CEP are canceled out, leading to a finite bulk viscosity 
at the critical end point [24(. 

The NJL model, being non-renormalizablc, requires 
the ultraviolet momentum cut-off A. Consequently, all 
thermodynamic quantities are affected by finite A. This 
is particularly the case in the parameter range where the 
physics is sensitive to a large particle momenta. Fig. [5] 
shows sensitivity of rj/s and £/s to the momentum cut- 
off At implemented in the thermal part of thermody- 
namic potential. We compare the results obtained with 
A T = A T=0 = 587.9 MeV and with A T =1 GcV. In 
high temperature phase we also show in Fig. [5] the re- 
sulting viscosities obtained from the perturbative QCD 
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0.7 0.8 0.9 1 1.1 1.2 1.3 0.7 0.8 0.9 1 1.1 1.2 1.3 

T/T c T/T c 

FIG. 8: The ratio of rj/s (left-hand figure) and £/s (right-hand figure) at fi = 0. The solid curves are obtained with the same 
3-momentum cutoff At in thermal part as in the vacuum, i.e., At = 587.9 MeV. The dashed lines are obtained for At = 1 GeV 
in themal part while in vacuum part A = 587.9 MeV is used. The dotted curves indicate the asymptotic values for massless 
quarks in high temperature limit. 



(pQCD) [H # 3 



and with A = 30 MeV. 



V 



(asymptotic) 



An 2 



3.86 



675 a 2 s ln(l/o! s ) 

^(asymptotic) q 



T 3 , 



with the two-loop running coupling for Nf = 2 



a* = 2n 



|ln(T/A) + ^ln(21n(T/A)) 



(3.15) 



(3.16) 



* 3 For the recent results of viscosi ty p arameters obtained in the 
high temperature QCD see Ref. l2- r l . 



As seen in Fig. [8l both rj/s and £/s are very sensitive 
to the momentum cut-off. However, at low T this de- 
pendence is opposite in shear and bulk viscosity. There 
is an increase of rj/s and decrease of £/s with increasing 
A. This is due to derivative terms in the bulk viscosity 
which show steeper change with A. At high temperature 
the NJL model results tend to converge to that expected 
in pQCD. The bulk viscosity vanishes at high tempera- 
ture, as expected in the conformal theory. However, due 
to the absence of gluons and the presence of the finite 
cut-off the shear viscosity deviates from their asymptotic 
value even at T ^> T c . 
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FIG. 9: The shear (left-hand figure) and bulk (right-hand figure) viscosities normalized by the 3-momentum cutoff calculated 
along the phase boundary. The vertical dotted-line indicates the position of the CEP. 
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FIG. 10: The ratio of shear (left-hand figure) and bulk (right-hand figure) viscosities to entropy density calculated along the 
phase boundary. The vertical dotted-line indicates the position of the CEP. 



3.2. Viscosities along the phase boundary 

Fig. [9] shows the shear and bulk viscosities calculated 
along the chiral phase boundary obtained in the NJL 
model and shown in Fig. [TJ Two lines at high /i c indi- 
cate a gap associated with the first order transition. The 
shear viscosity increases along the cross-over line towards 
the CEP where it reaches its maximal value. The bulk 
viscosity is rather weekly changing with T in the regime 
of cross-over and exhibits a local minimum at the CEP. 
A little bump on the cross-over line before reaching the 
CEP is an artifact of the present NJL model calculation. 
In a more realistic calculation, this may disappear and 
the bulk viscosity smoothly goes down to the CEP. The r\ 
shows an upward jump approaching the first order transi- 
tion from the low-temperature phase. This is an opposite 
behavior to £ which shows a downward jump. 



Fig. ITD1 shows the ratios, r\j s and £/s, calculated along 
the critical line. The r\j s is monotonic along the critical 
line without showing any sensitivity to the CEP. The 
C/s preserves the local minimum at the CEP and shows 
a visible gap at the first order transition. 

The properties of the transport coefficients calculated 
within NJL model are not universal in spite of the fact 
that this model belongs to the QCD universality class 
with respect to chiral symmetry In addition the shear 
and bulk viscosities are sensitive to particle composition 
of the medium and to specific model dynamics that enters 
to various collision processes. The shear and bulk vis- 
cosities do not show critical behavior with static critical 
exponents. Nevertheless, one expects that some generic 
futures of transport coefficients along the phase boundary 
calculated within the above model could be of interest for 
understanding transport properties of the QCD matter. 
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From this perspective, a trace of in-medium change of 
the order parameter M(T, /x) remains more distinctly in 
the bulk than in the shear viscosity. 

4. SUMMARY AND CONCLUSIONS 

Based on the NJL model and kinetic theory we have 
shown the transport properties of thermal medium com- 
posed of dynamical quarks. We have included in our 
calculations in-medium modification of constituent quark 
dispersion relations. We have analyzed behaviors of shear 
r] and bulk £ viscosities near the chiral phase transition 
within the mean field approximation. The collision time 
required to quantify viscosity parameters was calculated 
from the elastic scattering cross-section of constituent 
quarks including next to leading l/N c corrections. 

We have shown that around the phase transition the 
shear viscosity exhibits only a shallow minimum for 
small chemical potentials. For higher fi corresponding 
to the first order transition and near the critical end 
point the shear viscosity is discontinuous or exhibits an 
abrupt change. However, when normalizing rj by the en- 
tropy density s this abrupt change at CEP disappears. 
The bulk viscosity is more sensitive to the change of 
the dynamical quark mass around the phase transition, 
thus shows steeper behavior near T c , however stays non- 



singular at CEP. 

Along the chiral phase boundary the r\j s is changing 
monotonically with critical parameters whereas the Q/s 
has a local minimum at the CEP. Nevertheless, within 
the considered model and under relaxation time approx- 
imations none of viscosities could be considered as a use- 
ful probe of the critical end point in the phase diagram. 
We have to stress that due to the non-universal behavior 
of bulk and shear viscosity within these calculations our 
quantitative results arc model dependent. However, we 
expect that qualitative behavior of rj/s and £/s near the 
chiral phase boundary derived here could be less sensitive 
to the specific model dynamics since important change of 
the order paratmer is certainly incorporated in our for- 
mula. 
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